Theory of phonon-mediated relaxation in doped quantum dot molecules 
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A quantum dot molecule doped with a single electron in the presence of diagonal and off-diagonal 
carrier-phonon couplings is studied by means of a non-perturbative quantum kinetic theory. The 
interaction with acoustic phonons by deformation potential and piezoelectric coupling is taken into 
account. We show that the phonon-mediated relaxation is fast on a picosecond timescale and is 
dominated by the usually neglected off-diagonal coupling to the lattice degrees of freedom leading 
to phonon-assisted electron tunneling. We show that in the parameter regime of current electrical 
and optical experiments, the microscopic non-Markovian theory has to be employed. 
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I. INTRODUCTION 

A promising way of building feasible quantum com- 
puter is the control of quantum superpositions in cou- 
pled quantum dots (QDs) called quantum dot molecules 
(QDMs) 20 . In particular, the two electron ground states 
are used as the logical qubit states and can be manip- 
ulated by electrical 2 ^ as well as by optical mean o 22 ' 23 . 
However, the coherent control of the qubit states is lim- 
ited not only by the electron-electron interaction 2 ^ but 
also by the coupling to the semiconductor host materials 
resulting in emission and absorption of phonon o 25 ' 26 . 

Since the energy difference between the two electron 
ground states is of order of a few meV, this can lead 
to new phonon-mediated relaxation processes 2 ^. The 
electron transfer between the QDs is of great impor- 
tance for the currently performed optical and trans- 



port experiment 



,28,29 



but only little theoretical work 



exists on the influence of the phonon coupling. It 
was mostly limited to perturbative and Markovian 
theory or by including only diagonal electron-phonon 
interactio n 25 ! 30 " — . A comprehensive theory including 
non-Markovian processes and non-diagonal coupling to 
the phonon reservoir is required to adequately describe 
this important solid state qubit system. 

In this paper, we present a fully quantum kinetic de- 
scription of the phonon-mediated relaxation in doped 
quantum dot molecules including non-Markovian effects. 
We show that the phonon-electron interaction in cou- 
pled quantum dots leads to a fast electron relaxation on 
a picosecond timescale and strongly affects the coher- 
ent electron evolution. A non-perturbative theory based 
on correlation expansion techniqu o 33 ' including up to 
three particle correlations is used, which describes up 
to two-phonon assisted processes. We include both di- 
agonal and off-diagonal electron-phonon couplings and 
show that the usually neglected off-diagonal contribution 
can play a dominant role in relaxation processes leading 
to fast phonon-assisted electron tunneling. We derive a 
compact expression for the electron occupation and co- 
herence evolution in the Markovian approximation and 
compare the results with the full quantum kinetic de- 



scription. 

The paper is organized as follows. In Sec. [TTJ a quan- 
tum kinetic model describing an electron confined in a 
quantum dot molecule coupled to the phonon environ- 
ment is introduced. Next, Sec. IIIII contains the results 
on phonon-assisted electron evolution and the compari- 
son between non-Markovian and Markovian regime. In 
Sec. IIV1 we conclude the paper with final remarks. In 
addition, some technical details on the correlation ex- 
pansion technique are presented in the Appendix. 



II. MODEL SYSTEM 

The system under study consists of a singly doped 
quantum dot molecule build from two quantum dots. 
We consider the relaxation between the two energetically 
lowest states of the electron |1) in the left and |2) in the 
right quantum dot. The free electron Hamiltonian reads 

iJ c = e(|2)(2|-|lXl|)+r(|l}(2| + |2)(l|), (1) 

where T is the tunneling coupling between the QDs and 
2e is the energy difference between the ground states in 
both quantum dots. The free phonon Hamiltonian reads 
H p h — J2k s ^fc.s^fc s bk.s with phonon creation and an- 
nihilation operators b\. and bk, s , respectively, phonon 
frequencies uik.s and phonon wave vector k, where s la- 
bels different phonon branches (longitudinal I and two 
transversal tl and t2). The carrier interacts with the 
phonon reservoir, which is described by the following 
Hamiltonian 



Hint = 



E E 

i,J=l,2 k,s 
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where gij, s {k) are coupling elements of the electron- 
phonon interaction. Note, that we take into account the 
usually neglected off-diagonal electron-phonon coupling, 
describing the direct phonon-mediated electron transfer 
between the quantum dots with simultaneous emission 
or absorption of a phonon, which, as we will show, leads 
to a fast relaxation channel. 
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We consider the relevant acoustic phonons coupled via 
both piezoelectric coupling and deformation potential 
with corresponding coupling elements: 

&ik) = -i^^MM J *r (r)e^<,), 

(3) 



&) =\l^Dej d 3 r^(r)e ifer ^(r) (4) 
with Gaussian electron wave functions 

In the numerical simulations, the parameters correspond- 
ing to self-assembled GaAs quantum dots have been 
used: l\ = 4 nm, li = 4.1 nm, and d — 6 nm, 
where U is the electron wave function size of the i-th 
QD, and the distance between the dots is d. Here, 
p = 5360 kg/m 3 is the crystal density, V is the nor- 
malization volume of the phonon modes, c s is the speed 
of sound (longitudinal ci — 5150 m/s or transverse 
c t = 2800 m/s, depending on the phonon branch), dp = 
0.16 C/m 2 is the piezoelectric constant, D e = —8 eV 
is the deformation potential constant for the electrons 
and e r — 13.2 is the static dielectric constant. The 
function M s (k) in the piezoelectric coupling element de- 
pends only on the orientation of the phonon wave vec- 
tor and for a zinc-blende structure reads: Mi(k) = 
§sm(0)siri(20)sin(2v3), M n (fc) = - sin(20) cos(2ip), and 
M n (k) = sin(0)[3cos 2 (^) - 1] sin(2<p). 

In order to study the dynamics of the electron confined 
in a quantum dot molecule in the presence of electron- 
phonon interaction, a second order correlation expan- 
sion method is use d 33 ' 34 . This is a non-perturbative 
technique which can be employed for the microscopic 
quantum kinetic description of the dynamics covering 
the non-Markovian memory effects (see Appendix for 
details). The equations of motion for the electron oc- 
cupations (|2)(2|), coherences (|1)(2|), and the 
phonon-assisted elements e.g. (|l}{2|& (Jj , s ) co '" r are derived 
and solved by standard numerical methods. We included 
up to three particle correlations, e.g. (|1)(2|&£ s bk, s >) corr , 
which leads to a description of up to two-phonon assisted 
processes. An alternative approach to study the phonon- 
assisted dynamics is the real time renormalization-group 
method 3 ^. 



III. PHONON-ASSISTED ELECTRON 
EVOLUTION 

The electron is initially in a quantum dot with higher 
energy, / = (|2)(2|) = 1. We assume that the injection 
time is much shorter than the phonon reservoir response, 



so that all the correlation terms can be initially set to 
zero. In Fig. [TJ the occupation of the QD with higher 
energy with and without coupling to the phonon reservoir 
at four different temperatures is shown as a function of 
time. As expected, the ideal electron evolution without 
interaction with phonons (solid black line) is a coherent 
oscillation of the electron between the two quantum dots 
with period determined by the tunneling coupling and 
the energy difference between the two electron states Ae 
and is proportional to T/(Ae). For the presented set of 
parameters with Ae = 0.5 meV and T = 0.1 meV, the 
electron oscillates slightly but mostly stays in the QD 
with higher energy. The tunneling coupling is to small 
to let the electron fully tunnel between the dots. 

If the coupling to the phonon environment is included, 
the dynamics of the confined electron changes com- 
pletely: the electron tunnels to the neighboring dot al- 
ready on the picosecond time scale. At low tempera- 
ture, the electron tunnels from the quantum dot with 
the higher energy to the one with the lower energy. At 
T = K, only the emission processes are available, thus, 
the electron tunnels completely to the lower QD. As the 
temperature grows, the probability of the phonon ab- 
sorption increases. In consequence, the final state of the 
electron moves from the quantum dot with the lower 
energy towards an equal superposition of electron con- 
fined in the left |1) and right \2) quantum dot. This is 
due to the fact, that at high temperatures, the probabil- 
ity of phonon absorption equals to the probability of its 
emission. In addition, the amplitude of the oscillations 
decreases with growing temperature, since the damping 
gets stronger. Similarly, the initial decay is also faster at 
higher temperature. We compared the results including 
up to two-phonon assisted processes to those with up to 
one phonon and found out that the former leads only to 




time (ps) 



FIG. 1. The occupation probability of the quantum dot with 
the higher energy as a function of time at different temper- 
atures with Ae = 0.5 meV, and T = 0.1 meV, and the ideal 
electron evolution without coupling to phonons. 
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a small correction. Thus, the quantum kinetic descrip- 
tion with the first order correlation expansion technique 
results in the proper description of the electron dynamics. 



In order to show that the fast relaxation processes 
result from the direct phonon-assisted electron tunnel- 
ing described by the off-diagonal electron coupling to 
phonons, the occupation of a QD modeled using only the 
diagonal interaction Hamiltonian is shown in the inset in 
Fig. [2] One can see, that in this case the relaxation is 
considerably slower up to two orders of magnitude. Thus, 
in order to properly describe the dynamics of an electron 
in a quantum dot molecule, one cannot neglect the off- 
diagonal electron-phonon coupling. A similar conclusion 
has been met by Zibik et aL 36 ' 37 , who studied the co- 
herence of the s-p transition in single quantum dots. It 
has been shown that both direct absorption of acoustical 
phonons between the two p states as well as the scatter- 
ing via virtual transitions have to be included in order 
to explain the decoherence mechanisms of the considered 
transition. 



Next, for comparison with our full quantum kinetic 
description, we have derived a Markovian approximation 
for the electron dynamics governed by real electron tran- 
sitions, where the memory effects are neglected. Includ- 
ing phonon-assisted correlations up to the first order, it 
is possible to derive compact expressions for the electron 
density and coherence [see Eqs. © and ([7]), respectively]. 



The first term in the equation of motion for the density 
/ describes the ideal coherent evolution of the electron, 
whereas the second term corresponds to the correction 
due to the electron-phonon coupling. Here, one can eas- 
ily recognize the two terms corresponding to the phonon 
emission ~ f{n q , s + 1) and absorption (1 — f)n q , s pro- 
cesses, where n 9iS is the Bose distribution function. The 
occupation is directly affected only by the off-diagonal 
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FIG. 2. As in Fig. [T] but in the Markovian approximation. 
Inset: The evolution of the electron occupation including only 
the diagonal electron-phonon coupling. The lines and axes as 
in Fig. [Lj 
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FIG. 3. The time evolution of the electron density in the QD 
with higher energy for different energy differences between the 
electron states with T = 0.1 meV and at T = 4 K. 



coupling to phonon reservoir. In contrast to this, the 
polarization contains also the diagonal elements which 
lead to pure dephasing processes 3 ^—. They are particu- 
larly visible in the case of quantum dots with large size 
difference. It is worth mentioning that the considered 
phonon-induced decay is a few orders of magnitude faster 
than the radiative decay of an exciton confined in a sin- 
gle quantum dot, which for self- assembled GaAs QDs is 
in the order of 1 ns^i. The results of the numerical sim- 
ulations for the Markovian approximation are shown in 
Fig. [3 Here, the coherent oscillations of the electron 
vanish due to strong diagonal and off-diagonal damping 
elements. Thus, the difference between the quantum ki- 
netic description and Markovian approximation are par- 
ticularly visible in the case of lower temperatures, where 
the amplitude of the oscillations is large in the full de- 
scription. 

As mentioned previously, the tunneling coupling and 
energy difference determine the period and amplitude of 
the coherent oscillations between the two electron states 
in both quantum dots. In Fig. [3l the time evolution of 
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./ = |rim(p) + l5i2,,(g)| 2 5 (fiw,, a - 2e) {2 [/(n g , s + 1) - (1 - /)»»,,.] - 4Im 2 (p)} 



(6) 



P = -r2ep - ~r(l - 2/) + J 5> 2 \ s (<z) - <7n, s (<z)]Si2 >s OK,. - 2e) [/(n,, s + 1) - (1 - f)n q , 8 - \ P \ 2 ] (7) 



q,s 



-\ ^2[922, s (q) ~ 5ii,s(g)]5i*2, s (9)p 2( ^ {f^q,s - 2e) - | ^ I5i2,s(g)| 2 (/ + ng :S )4iim(p)<S (fiw gjS - 2e) . 
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FIG. 4. As in Fig. [3] but with the Markovian approxima- 
tion. Inset: The ideal evolution without interaction with the 
phonon reservoir. 



If the Markovian approximation is employed [see 
Fig. 2] the coherent oscillations vanish, thus the discrep- 
ancies from the quantum kinetic theory is largest for the 
energy differences approaching the value of the tunneling 
coupling, (the ideal electron evolution is shown in the in- 
set in Fig. 2] for reference). This case corresponds to the 
regime which is common in the recently performed trans- 
port and optical experiments. Thus, in order to properly 
describe the electron dynamics, the full non-Markovian 
description including the off-diagonal coupling to the lat- 
tice degrees of freedom has to be used. In the case of 
energy differences much larger compared to the tunnel- 
ing coupling, the Markovian approximation delivers the 
same results as the quantum kinetic description. 



IV. CONCLUSION 

In conclusion, we have presented a full description of 
the electron dynamics in the presence of the electron- 
phonon coupling in a quantum dot molecule doped with 
a single electron. This solid state system is of great 
importance for the nanoscience, in particular for quan- 
tum information processing as well as quantum transport 
and optics. The two couplings to the relevant acoustic 
phonons via deformation potential and piezoelectric cou- 



I 

pling have been taken into account. We showed that 
the phonon-mediated relaxation is a fast process on a 
picosecond timescale and strongly modifies the coherent 
evolution of the electron. Based on the full quantum ki- 
netic description, we derived a compact expressions for 
the electron density and polarization in the Markovian 
approximation and we showed that in the regime of the 
current electrical and optical experiments, the full non- 
Markovian description is needed. We showed that the 
usually neglected off-diagonal electron-phonon coupling 
is dominant in the electron dynamics and leads to two 
orders of magnitude faster relaxation. 



ACKNOWLEDGMENTS 

A. G.-G. thanks P. Machnikowski and C. Emary for 
helpful discussions. A. G.-G. and J. F. acknowledge 
support from the Emmy Noether Program of the DFG 
(Grant No. FO 637/1-1) and the DFG Research Train- 
ing Group GRK 1464, and thank the John von Neumann 
Institut fur Computing (NIC) for computing time. 



Appendix: Correlation expansion technique 

In this Appendix, we present some details concerning 
the correlation expansion techniqu e) 33 ' 34 , which is based 
on the assumption that correlations involving an increas- 
ing number of particles are of decreasing importance. It is 
a non-perturbative method widely used for microscopic 
quantum kinetic description of the dynamics including 
the memory effects in the non-Markovian regime. 

Using Heisenberg equation one derives the equations of 
motion for the quantities of interest: electron occupation 
/ = (|2)(2|) and coherence p = (|1)(2|). The equation for 
the time evolution of the electron density 



/ = jT Im(p) - | 



E^( fc ) 



Sfe 



„(+)* 



2iB k lm(p) 



(+) 



2iB* k lm(p) 



couples to the coherence p, B k — (bk), and to the pho non- 



assisted correlations Sk 



(\m\bkY 



and 



„(+) 



(|l)(2|6l) c< 

used: (|l)(2|h fc ) = (|l)(2|)(6 fe ) + (|1}<2|6* 



Here, the factorization scheme has been 

r , where the 
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quantities have been decompose into all possible lower- 
order factorizations. The next step is to derive the equa- 
tions of motion for p and all the phonon-assisted cor- 
relation quantities. The latter couple to three particle 
correlations, e.g. (|1)(2|&1 s bk :S ') corr which couple to four 



particle correlations, etc. Thus, in order to get a closed 
set of equations, one needs to truncate the hierarchy by 
neglecting higher order correlations. In the present pa- 
per, we included up to three particle correlations. A com- 
prehensive review on the correlation expansion technique 
is presented in Ref. HH- 
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